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ABSTRACT 


A numerical study was conducted to predict the 
combined convective and radiative heat transfer rates on 
the walls of a small aspect ratio cylinder representative 
of the scaled model of a rocket engine combustion chamber. 
A high-temperature, high-pressure environment was simulated 
in the cylinder, with gas velocities at low subsonic levels 
typical of the conditions leading to the entrance of the 
nozzle section of a rocket engine. The composition of the 
gases in the cylinder was determined from the TEP program 
for the burning of rocket fuel at typical values of the 0/F 
ratio. The thrust of the study was to determine the 
radiative contribution to the heat transfer rate from the 
hot participating chamber gases to the cooler wall. The 
calculations were carried out using the commercial CFD 
package CFDACE, and were first benchmarked against known 
results in the literature for the simpler case of gray 
chamber walls and a gray participating medium. The non-gray 
computations were subsequently carried out using gas 
absorption coefficient values obtained from the exponential 
wide band model with the help of the fire-modeling program, 
RADCAL. The effect of different chamber wall temperatures 
and gas compositions were examined. The main findings of 
the study are that the radiative contributions at the high 
gas temperatures being considered are comparable to the 
convective values, and strongly spectral in nature. In 
addition these radiative fluxes were found to be least 
sensitive to the wall temperature and chamber pressure in 
the range considered. Furthermore, this radiative 
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contribution reaches a maximum at a unique optimal optical 
thickness of the gas that lies within the extremes of the 
optically thin and thick limiting cases. 
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I. 


INTRODUCTION 


The heat transfer behavior in many high temperature 
applications involves a combination of radiative and 
convective heat transfer modes. Analyses of such 
applications are often complicated by irregular geometry, 
turbulence, and combustion. 

Radiation within a medium containing products of 
combustion is dependent upon the temperature and 
concentrations through the entire field. The energy is 
distributed across the infrared spectrum in a highly 
nonlinear fashion. This greatly complicates modeling of the 
heat transfer within such high temperature environments. 

Gaseous radiation is particularly relevant to heat 
transfer in combustion chambers burning carbon, hydrogen, 
or hydrocarbon fuels. Absorption and emission by carbon 
dioxide, water vapor and other products of combustion is of 
particular concern due to their strong radiative 
participation. 

This thesis deals with the characteristics of combined 
convective and radiative heat transfer rates on the walls 
of a small aspect ratio cylinder representative of the 
scaled model of a rocket engine combustion chamber. 
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II. BACKGROUND & MOTIVATION 


The flow and heat transfer behavior in a rocket engine 
is crucial in determining its performance as well as its 
thermal signature. 


The schematic of a typical rocket engine is shown in 
Figure 1. 



Figure 1. 


Rocket Engine Schematic 


There are three main sections in a rocket engine. The 
first is the combustion chamber, in which oxygen as 
oxidizer and fuel, i.e. RP1, are ignited and a high- 
temperature, high pressure environment is created. Gas 
velocities are at low subsonic levels typical of the 
conditions leading to the nozzle section. In the nozzle the 
flow is accelerated to supersonic velocities and is ejected 
in the exhaust plume. 


For 

this 

phase of 

the 

research program. 

this thesis 

focuses 

only 

on the 

heat 

transfer behavior in 

the 

combustion chamber and 

this 

thesis provides 

modeling 

and 
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III. PROBLEM DESCRIPTION 


A. THE MODEL PROBLEM 

Detailed schematic of scaled-down model of cylindrical 
combustion chamber as used in the experimental program is 
shown in Figure 2. 
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The scaled model of the combustion chamber consists of 


a participating (absorbing and emitting but non-scattering) 
gray/non-gray gas flowing through a circular duct in the 
entrance developing region of flow. The gas is assumed to 
enter the duct at a uniform velocity and high temperature 
and pressure. The duct wall is treated as isothermal with 
temperatures ranging from 400 °K to 800 °K . A schematic 
diagram of the model is shown in Figure 2.1. The inner 
diameter of the duct is 1.5*, and its length is 3*. RP1 is 
the fuel being burned, while oxygen is used as oxidizer. 
The products of the combustion process are determined based 
on an O/F value of 3 which is typical of the chamber. The 
actual combustion is not modeled in this initial study, 
although its net effect is represented as a well mixed core 
of the constituents with core temperatures ranging from 
1500°K to 3500°K for use in this parametric study. The 
inlet fluid velocity to the duct is taken to be uniform at 
a value corresponding to a low subsonic Mach number. 

B. GAS RADIATION 

As per the laws of quantum mechanics, molecular gases 
can emit or absorb photons to varying degrees at an 
infinite set of distinct wavenumbers or frequencies. No 
spectral line can be truly monochromatic; rather, 
absorption or emission occurs over a tiny but finite range 
of wavenumbers. The results are broadened spectral lines 
that have their centers (maxima) at the wavenumber 
predicted by quantum mechanics. Three most important 
phenomena that cause broadening of spectral lines are 
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collision broadening, natural line broadening, and Doppler 
broadening. 

A single spectral line at a spectral position is 
characterized by its intensity and its line half-width. 
However, a vibration-rotation band has many closely spaced 
spectral lines that may overlap considerably. The 
absorption coefficient for the entire band is found by 
adding the absorption coefficients for single lines. 

= IX- (2 - 1} 

j 

All spectral integrations may be reduced to two cases, 

1 K vhb)r, dT l and | 0 Cw[ | -exp(-{V, ; ds')]J 7 (2.2) 


where 7 (i) denotes that either I b or I can occur. 


Since the local radiation intensity is due to emission 
from all locations within the medium and the bounding 
walls, it can be expected to be smooth. It is further 
smoothed by absorption and scattering. As a result, the 
previous equation simplifies to 



1 M7+A17/2 

[At] 2 


K,dV 


dr i 1 


(2.3) 


and 




(2.4) 


In order to find narrow band values of the absorption 
coefficient and the emissivity, required information about 
the spacing of the individual lines within the group and 
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their relative strength is supplemented by a number of 
models. Two extreme models are Elsasser model, in which 
equally spaced lines of equal intensity are considered and 
statistical model, in which the spectral lines are assumed 
to have random spacing and/or intensity. A typical spectral 
line arrangement for these models is shown in Figure 3. 



Figure 3. Typical Spectral Line Arrangements for 

(a) Elsasser and (b) Statistical Model [1] 


In statistical model it is assumed that the spectral 
lines are not equally spaced but, rather, are randomly 
distributed across the narrow band. A statistical analysis 
can be carried out in two different ways. In the uniform 
statistical model the line intensity is assumed to be 
constant, leading to 


W 

£ = 1 - exp(—) 
d 


where 


W 

d 


(2.5) 

-spectrally averaged emissivity 

-equivalent line width 
-spacing between spectral lines 



In general statistical model, the statistical model 
has also been evaluated for exponentially decaying line 
intensities, resulting in 


= 1 - CX P (_ "T==) = 1 “ ex P( ) 

i T 

1 +- 1 + — 


( 2 . 6 ) 


where: £ -spectrally averaged emissivity 


S -line-integrated absorption coefficient 

d -spacing between spectral lines 

X -optical path length 

b L -line half width 


-line overlap parameter 

-narrow band optical thickness 

The results from statistical model is summarized in 
Figure 4 and Table 1 in [1] as 
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Weak line Strong line No overlap All regimes 


x«l 


line, W SX 


x > 1 


0«1 


Elsasser Model, e n 
S - const 1 -1 
d - const 


Statistical Model, e ^ 

S = const 1 - e 
d random 





S exponential 1 - e 7 1 - expf- 


SXe x [k M + /|(x)] 


2 T A[\- e ~**M) 


er T7y 


/ W 

l-ex P - 7 


d random 


Definitions: x ~ fi = ir—, r = -X - 2 fix 
2vb L d d 
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Line overlap parameter /3 



0.001 


0.0001 - 

-1-1__lJ_I_ 

0.001 0.01 0.1 1.0 10 100 

Line strength parameter x 

Figure 4. Regions of Validity for Elsasser and 

Statistical Band Model Regimes [1] 


C. RADIATIVE TRANSPORT EQUATION (RTE) 

In evaluating the radiative energy transport, it is 
required to know radiative properties, such as, emissivity. 


11 



absorptivity, and reflectivity in the case of surfaces, and 
absorption coefficient and scattering coefficient for the 
participating media such as gases. The law of conservation 
of energy is then applied to determine the energy field. 

The general relationship that governs the behavior of 
radiative heat transfer in the presence of an absorbing, 

emitting medium is developed in [1] . RTE describes the 
radiative intensity field within the enclosure as a 

function of location, fixed by location vector r, 

direction, fixed by unit direction vector , and spectral 
variable, wavenumber . The net radiative heat flux 

crossing a surface element is obtained by adding up the 

contributions of radiative energy irradiating the surface 
from all possible directions and for all possible 
wavenumbers. Integrating the equation of transfer over all 
directions and wavenumbers leads to a conversation of 
radiative energy statement applied to an infinitesimal 

volume. 

If the medium is participating, in that case any 
incident beam will be attenuated by absorption and 
scattering while radiative energy travels the medium, as 
shown in Figure 5. 

The amount of absorption is directly proportional to 
the magnitude of the incident energy as well as the 
distance the beam travels through the medium. Thus, 

( dI V )ab S =- K V I V ds ( 2 - 7 ) 
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is the linear 


where the proportionality constant (7 sn 

scattering coefficient for scattering from the pencil of 
rays under consideration into all other directions. 

A light beam traveling through a participating medium 
in the direction of s loses energy by absorption and 
scattering away from the direction of travel. But, it also 
gains energy by emission and by scattering from other 
directions into the direction s . 

The rate of emission from a volume element is 
proportional to the magnitude of the volume. The emitted 
intensity along any path is proportional to the length of 
the path and the local energy content in the medium. At 
thermodynamic equilibrium, the intensity everywhere must be 
equal to the blackbody intensity. 

(dl n ) em = K n I bJ1 ds (2.9) 

Augmentation due to in-scattering has contributions 
from all directions, and it is calculated by integration 
over all solid angles. The energy flux coming from all 
directions and scattered into the direction s is given in 
[1] as 

( dI n ) sca =ds °^j. L ! n ( S n S ) dSl (2.10) 

An energy balance from equations (2.7), (2.8), (2.9), 

and (2.10) results to RTE in [1] as 

I l/ (s+ds,sj+ dt ) -I n (s, s, t ) = K n I hn (s, t)ds - K n l q (s, s, t)ds - (T slJ I l} (s, s, t)ds 

a , ( 2 . 11 ) 
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D. 


THE METHOD OF DISCRETE ORDINATES (S n -APPROXIMATION) 


The Discrete Ordinance Method (DOM) is a tool to 
transform the equation of transfer into a set of 
simultaneous partial differential equations. DOM is based 
on a discrete representation of the directional variation 
of the radiative intensity. A solution to the transport 
problem is found by solving the equation of transfer for a 
set of discrete directions spanning the total solid angle 
range of 4n . The discrete ordinates method is simply a 
finite differencing of the directional dependence of the 
equation of transfer. Integrals over solid angle are 
approximated by numerical quadrature. 


The general equation of transfer for an absorbing, 
emitting, and anisotropically scattering medium is 


dl 

ds 


= s VI(r,s) = fc(r)I b (r)~J3(r)l(r,s) + 


4 n 


J I(r, s', s)dQ.' 

4 Jt 


(2.7) 


Equation (2.7) is valid for a gray medium or, on a 
spectral basis, for a nongray medium, and is subject to the 
boundary condition 

I(r w , s) = £{r w )I b (r ) + f.. I(r w ,s)\n-s\fin' (2.8) 

Jns< 0 1 1 

where the enclosure is considered as opaque, diffusely 
emitting and reflecting wall. 

In the discrete ordinates method, equation (2.7) is 
solved for a set of n different directions s jf i = l,2,0,n, 

and the integrals over direction are replaced by numerical 
quadratures, that is. 
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(2.9) 


£ f(s)dCl = 

n i=i 

where, w t is the quadrature weights associated with the 
directions S; . 

And, equation (2.7) becomes, 

(7 (r) " 

■ V/(r, s ,) = K{r)I b ( r ) - J3(r)l(r, s ,) + — ^ Wjl(r, s . )4>(r, s,, s .), (2.10) 

4 ^ Jl J 

where, j=l,2,0,n, subject to boundary conditions 

I ( r w ^ i ) = £{r w ) I b {r w ) + ^^- 'E i w J I(r w ,s J )\n-s J \ (2.11) 

71 «Sj< 0 

After determining the intensities, the radiative heat 
flux inside the medium or at a surface may be found as 

n 

q(r) = [ I(r, s)sd£l = Y w i I I (; r)s i (2.12) 

J47T ' 4 

i =1 

In Table 3.2, the direction cosines of s ; and 
corresponding weights for S 4 -approximation is shown. The 
direction cosines of s ; are 

^ = (s t ■ i )i + (s t ■ j)j + (s t ■ k)k = + rjJ + //,i (2.13) 

Only positive direction cosines are given in Table 
3.2, covering one eight of the total range of solid angles 
4 . To cover the entire 4 any or all of the values of^, ri i( 
and may be positive or negative. 

Only positive direction cosines are given in Table 2, 
covering one eight of the total range of solid angles 4 . 
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To 

and 


cover the entire 4 any or all of the values of^ ; , 
(. 1 ; may be positive or negative. 


Order of 
Approximation 

Ordinates 

Weights 


5 

n 

h 

w 

S 4 

0.2958759 

0.2958759 

0.9082483 

0.5235987 


0.2958759 

0.9082483 

0.2958759 

0.5235987 


0.9082483 

0.2958759 

0.2958759 

0.5235987 


Table 2. Discrete Ordinates for the S 4 


Approximation [1] 
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IV. NUMERICAL TOOLS 

A. NUMERICAL SOLUTION STRATEGY 

There are mainly three numerical tools used in this 
study. 

First tool is TEP, which is used to calculate the 
thermodynamic chemical equilibrium properties of gas 
mixtures resulting from the combustion of various fuels and 
oxidants. 

The second one is RADCAL, which is a numerical fire¬ 
modeling program which provides detailed gas spectral 
radiation data. 

The third one is CFDACE, which is a set of computer 
programs for multi-physics computational analysis. 

The numerical solution strategy flowchart is shown in 
Figure 6. 



Figure 6. Numerical Solution Strategy Flowchart 

B . TEP (THERMO-CHEMICAL EQUILIBRIUM PROGRAM) 

TEP is used to calculate the thermodynamic chemical 

equilibrium properties of gas mixtures for a variety of 

different applications. Both thermodynamic data and 

ingredient libraries are included with TEP. The 
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applications that are included in the TEP program are 
Rocket Application, Gas Properties, Shock Application, 
Scramjet Application, Isentropic Path, and Detonation 
Application. 

The application that is used is rocket problem. The 
rocket problem calculates the "theoretical", or maximum 
attainable performance, of a rocket thrust chamber. The 
calculations are based on the following assumptions: 

• One-dimensional form of the continuity, energy, and 
momentum equations 

• Zero velocity (stagnation) in the combustion chamber 

• Complete combustion 

• Adiabatic combustion 

• Isentropic expansion 

• Homogeneous mixing 

• Ideal gas law 

• Zero temperature and velocity lags between condensed 
and gaseous species 

For equilibrium performance, composition is assumed to 
reach equilibrium immediately during expansion. For frozen 
performance, composition is assumed to remain fixed at the 
initial (stagnation) combustion composition during 
expansion. 

Rocket performance is obtained by first determining 
stagnation combustion properties in the rocket chamber 
based on input of both enthalpy and pressure (H, P) or 
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temperature and pressure (T, P). Next, conditions at the 
nozzle throat are determined. The procedure used here is 
to find the pressure at which the flow velocity is equal to 
the sound speed. Iteration is required using entropy, 
pressure (S, P) point solutions, where the entropy is held 
fixed at the value calculated previously in the chamber. 
Once the throat conditions have been determined, solutions 
at input area ratios are obtained by iteration using (S, P) 
point solutions. For a given area ratio, A/A*, the entropy 
is held constant, and pressure is iterated until the 
product of density and velocity satisfies the continuity 
relation; i.e., a value equal to the throat density x 
throat velocity/ (A/A*). 

There are two main inputs to the rocket application, 
thermochemistry and rocket operating conditions. In 
thermochemistry chemical reactants, thermochemical species 
data, and reactants mixture specification is defined. Under 
thermochemistry chemical reactants option, the ingredients 
that will be used as oxidizer and fuel are selected. And 
oxidizer to fuel ratio is entered. Under thermochemical 
species data option, the species that will be used in the 
application is selected between CHNO, CHO, HO, SRM, Master 
Thermo, and Air/Hydrogen file. The properties of these 
files are as follows: 

CHNO Thermodynamic Data File contains 29 species 
screened for use in air/hydrocarbon combustion problems. 
This set of species is more than sufficient for most CHNO 
problems. The study was done by screening calculated 

results obtained using the master files. The calculations 
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were for air/acetylene and air/methane for a broad range of 
pressure, temperature, and equivalence ratio. 

CHO Thermodynamic Data File contains 20 species 
screened for use in oxygen/hydrocarbon combustion problems. 
It has been obtained from the CHNO thermo data file by 
omitting species containing N. 

HO Thermodynamic Data File contains 8 species screened 
for use in oxygen/hydrogen combustion problems. It has 
been obtained from the CHO thermo data file by omitting 
species containing C. 

Solid Propellant Rocket Motor File contains the 
principal species that make up the exhaust of aluminized 
solid propellant rocket motors. 

Master Thermodynamic Data File contains data for 
approximately 1000 chemical species. It can be edited to 
produce short files containing only those species important 
to specific problems. Several files of this type are also 
provided. Principal species of the following elements are 


contained in 

this master 

f ile: 

Al, 

Ar, 

B, 

Ba, 

Be, 

Br, 

c. 

Ca, Cl, Cr, 

Cs, Cu, D, F, 

Fe, H, 

He, 

Hg 

, I, K, 

Li, 

Mg, 

N, 

Na, Nb, Ne, 

Ni, 0, P, Pb, 

S, Si, 

Sr, 

T a, 

Ti, 

V, 

Xe, 

Zn, 

and 


Zr. The electron and positive or negative ion forms for a 
limited number of species are also provided. Please note 
that the master file is intended as a library, and not for 
general problem solving use. The master file will often 
yield too many species for problem solving, especially when 
a reactant contains the element carbon. 
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Air/Hydrogen File is intended for use in Air/Hydrogen 
combustion problems, such as occur in flame propagation and 
shock tube studies. 

Under rocket operating conditions option, chamber 
pressure in atmosphere, chamber temperature in Kelvin, 
subsonic area, supersonic area, and pressure ratios are 
entered. 


Table 3 shows the input and output data for TEP 
program. 


INPUT DATA 

OUTPUT DATA 

Gas temperature, gas 

Chamber and exit 

pressure, oxidizer to fuel 

pressure, participating gas 

ratio, subsonic and 

temperature, density. 

supersonic area ratios 

molecular weight, sonic 

velocity, prandtl number, 

mixture species! mole 

fractions, mass fractions. 


Table 3. Input and Output Data for TEP Program 


C. RADCAL 

RADCAL is a narrow band numerical model developed at 

NIST for radiation calculations in a combustion 

environment. RADCAL is a numerical program, which predicts 

the radiant intensity leaving a nonisothermal volume 

containing nonuniform levels of carbon dioxide, water 

vapor, methane, carbon monoxide, nitrogen, oxygen, and 

soot. The absorption coefficient of the combined gases is 
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calculated from a narrow-band model and a combination of 


tabulated spectral properties and theoretical 
approximations to the vibrational-rotational molecular 
bands. 

Table 4 summarizes the species and their molecular 
bands currently in RADCAL. 


Species 

Band 

C0 2 

2.0 pm 

co 2 

2.7 pm 

co 2 

4.3 pm 

co 2 

10.0 pm 

co 2 

15.0 pm 

H 2 0 

1.38 pm 

H 2 0 

1.88 pm 

H 2 0 

2.7 pm 

H 2 0 

6.3 pm 

h 2 o 

20-200 pm 

CO 

4.6 pm 

ch 4 

2.4 pm 

ch 4 

3.3 pm 

ch 4 

7.7 pm 

soot 

0.4-2000 pm 


Table 4. Molecular Bands Included in RADCAL [3] 

The numerical code consists of a short main program 
which reads the temperature and concentration information 
from the data file, RC DAT. SUBROUTINE RADCAL is called 
from the main program to perform all the calculations. 
RCOUT DAT is the file that contains the outputs. Figure 7 
is a diagram of the structure of the program. 

RADCAL relies upon the four subroutines to compute the 
narrow-band parameters for the carbon dioxide 
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(SUBROUTINE C0 2 ) , water vapor (SUBROUTINE H 2 0 ) , carbon 
monoxide (SUBROUTINE CO ) , and methane (SUBROUTINE FUEL). Two 
large block data files (BLOCK DATA BD1 and BLOCK DATA BD2) 
contain the absorption coefficient of water vapor as a 
function of temperature and wave number. A third data file, 
BLOCK DATA BD3, contains similar information for the 15 m 
band C0 2 , and the 3.3 and 7.4 m bands of CH 4 . 

When the path contains only one element three 
different absorption coefficients are calculated the 
Planck-mean, the incident-mean, and the effective 
absorption coefficient. 

In the first line of the input data file, RC DAT, the 
number of elements into which the path is divided, NPT, is 
written. The second line lists the size of the first 
element in meters, DD (1), its temperature in Kelvin, T 
(1), and the partial pressures in kilopascals of carbon 
dioxide, P (1,1), water vapor, P (1,2), methane, P (1,3), 
carbon monoxide, P (1,4), oxygen, P (1,5),and nitrogen, P 
(1,6) . The last entry on the second line is the volume 
fraction of soot in the first element, W (1) . The third 
line contains the size, temperature, species partial 
pressure, and soot volume fraction for the second element. 

Similar data is entered for the remainder of the NPT 
elements. Following this information is a line containing 
the wall temperature in Kelvin, TWALL, and the minimum and 
maximum wavenumbers incm" 1 , OMMIN and OMMAX. The program 
terminates when it reads zero for the value of NPT. 
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Figure 7. Structure of Radiation Calculation 

Program RADCAL [3] 


The results of the calculations are found in the 
output file, RCOUT DAT. The input conditions are summarized 
in tabular form, followed by the total directional radiated 
energy flux emanating outward from element one, Q. The 
spectral intensity, QW (K) and the transmittance, TTAU (K) , 
are listed for each wavelength, AMBDA (K). The program also 
calculates the effective-mean absorption coefficient, 
AMEAN, the Planck-mean absorption coefficient, APO, and the 
wall-incident-mean absorption coefficient, AIWALL. 
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Table 5 shows the input and output data for RADCAL 


program. 


INPUT DATA 

OUTPUT DATA 

The 

number of 

elements, 

Total 

directional 

size i 

of the 

elements, 

radiated energy 

flux, the 

temperature of the 

elements, 

spectral intensity and the 

partial 

pressures 

of carbon 

transmittance 

for each 

dioxide, 

water 

vapor. 

wavelength, the 

effective- 

methane. 

carbon 

monoxide. 

mean absorption 

coefficient , 

oxygen. 

and nitrogen, volume 

the Planck-mean 

absorption 

fraction 

of soot, wall 

coefficient, and 

the wall- 

temperature, minimum and 

incident-mean 

absorption 

maximum 

wavenumbers 

• 

coefficient. 



Table 5. Input and Output Data for RADCAL 


Program 


D. CFD-ACE+ 

CFD-ACE+ is a set of computer programs for multi¬ 
physics computational analysis released by CFD Research 
Corporation. The programs provide an integrated geometry 
and grid generation module, a graphical user interface for 
preparation of the model, a computational solver for 
performing the simulation, and an interactive visualization 
program for examination and analysis of simulation results. 
CFD-ACE+ package include the following applications: 
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• CFD-GEOM Interactive 3D Geometry Modeling and 
Mesh Generation (Structured, Unstructured, and Mixed- 
Element Meshes) 

• CFD-GUI Solver Setup Interface 

• CFD-ACE (U) General Purpose Multi-Disciplinary 
Physics Solver 

• CFD-VIEW Interactive 3D Graphics, Animation & 
Flow Visualization Software 

Figure 8 shows the applications that are included in 
CFD-ACE + package. 



Figure 8. CFD ACE+ Package [9] 

1 . CFD-GEOM 

CFD-GEOM is an interactive CAD type geometry creation 
and fast grid generation (structured, unstructured and 
hybrid grids) program for Computational Fluid Dynamics 
(CFD). The geometry module allows the user to interactively 
define the geometry or modify the geometry read through the 
Plot3D data format. 

It serves as the front-end of CFDRC's solver packages 


and can be used with any other CFD/FEA/CAE program. With 
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interfaces to read geometry data created in almost any 
major commercial CAD program, it integrates CFDRC's CFD 
packages into the existing CAE environment. 

CFD-GEOM uses efficient data structures providing 
flexible inter-domain connectivity, part merging, and easy 
manipulation of all entities (geometry, topology, 
structured, unstructured meshes) with rapid database 
updating. Entity manipulations include translation, 
rotation, scaling, copying, (un)deleting, blanking, 
labeling, and model merging/extracting. Topological 
entities include: 

• Edges - composed of single or multiple curves and/or 
lines. 

• Faces - comprising single or multiple edge sets or 
single surfaces, possibly projected to surfaces. 

• Blocks - consisting of single or multiple volumes or 
face sets. Arbitrary orientation between adjacent 
blocks is allowed, or a consistent multi-block 
orientation can be requested automatically (hence no 
need to pre-plan individual IJK node directions). 

• Composition of edges, faces and blocks. 

• Edge-linking - any number of edges can be linked to a 
smaller number of master edges for mesh density and/or 
distribution, requiring minimum effort to modify mesh 
characteristics. 

• Boundary Conditions- the location and type of boundary 
condition for the flow solvers can be prescribed in 
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CFD-GEOM, providing automatic updating during mesh 

refinement or geometry change. 

2 . CFD-ACE (U) 

CFD-ACE is an unstructured, polyhedral cell flow 
solver. It is also integrated with a wide variety of other 
physics modules making it the core of a multi-disciplinary 
analysis environment. CFD-ACE employs a cell-centered 
control volume solution approach. This approach implies 
that the discrete equations are formulated by evaluating 
and integrating the fluxes across the faces that surround 
each control volume. 

CFD-ACE provides an integrated geometry and grid 
generation module. 

3 . CFD-VIEW 

CFD-VIEW is a 3-D post processor included in CFD-ACE+, 
which reduces the large volume of data generated by CFD-ACE 
to useful information. CFD-VIEW is an interactive graphics 
program for post-processing numerical results from CFD and 


other analysis 

software. 

It 

provides 

an easy- 

-to-use 

and 

interactive environment. 

with many 

graphics 

tools 

to 

visualize the 

flow physics. 

as well 

as the 

ability 

to 


extract data relevant to engineering design. 
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V. 


CODE CALIBRATION AND BENCHMARKING 


A. NUMERICAL TOOL CALIBRATION 

Geometry of the model is created by the pre-processor 
CFD-GEOM. Creating the structured grid is done by grid 
generation module in a bottom-up approach in labeling 
geometry elements as edges, faces, and blocks. In order to 
achieve the best results, the calibration is done when the 
participating medium is not present. This helps to compare 
the results with the hand calculations. 

The following parameters are considered when 
calibrating the model: 

• The edge parameters 

• The number of iterations 


• The error in the wall convection 


The model is created by 25, 50, 100, and 200 points 
consecutively. The edge parameters dictate how many grid 
points are placed on the edge and how these grid points are 
distributed. The power law is used as a grid point 
distribution method. The power law distribution is defined 
by 


u(n) = 


f 

V 


n — 1 

npts -1 


Y 

) 


(5.1) 


where: u(n) -location of grid point on the 

interval [0 01] 
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n -grid points index on the interval 

[ lOnpts] 

npts -the number of grid points 


x -user-specified power 

The power law is usually close to one which 
corresponds to a uniform distribution. For non-unity 
powers, the forward, backward, and symmetry options 
indicate how the power law grid distribution is applied. 


The geometry is created with 25, 50, 100 and 200 

points, while the power law grid scaling scheme is varied 
between 1.3, 1.5 and 2. Convergence criteria were varied 

from le-04 to le-10. Figure 9 discusses the results for 
convergence criteria of le-08. 
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Figure 9. Total Heat Flux Values for Various 

Power Law Grid Scaling Schemes 

As seen in the figure, the heat flux value for 100 
points with power law of 2 converges more rapidly than the 
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others and a 100 point grid was found to be a good trade¬ 
off between accuracy and computational time. 

B. TRASPARENT MEDIUM CASE 

Numerical results are compared with the analysis of 
radiative heat transfer within enclosures without 
participating medium. Figure 10 defines the surfaces and 
the emissivities corresponding to these surfaces. 


p T 



Figure 10. Two Dimensional Cylinder 

Radiative exchange between black surfaces is 
calculated from the following expression in [1] as 

q i =f J F i -AE hi -E bj ) (5.2) 

i=i 

where: 

N -number of surfaces 

F t _j -view factor from surface i to j 

E hi -blackbody emissive power of surface i 

The equation 5.2 becomes 
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‘h — ^ 1-1 (E bX E hX ) +F X _ 2 (E M E b2 ) +F x _ 3 (E bl E h3 ) 


(5.3) 


where the view factors can be calculated from view factor 
catalogue. [2] 


Table 6 shows the comparison of the hand calculations 
and CFD-A software results for wall temperature of 400 F 
and gas temperature of 3500 K between black surfaces. 



Hand 

Calculation 

Numerical 

Results 

Relative Error 

(%) 

Radiation to the 




Wall (watts/rad) 

-2882.46 

- 

3.02 



2795 



Table 6. Relative Error in Radiative Exchange 

Between Black Surfaces 


Similarly, radiation exchange between gray and diffuse 
surfaces can be calculated from the following expression 
[ 1 ] 


c 1l 


N 1 

K-- 1 ^ 






Z F i-j( E »- £ b j) 


7=1 


(5.4) 


The equation (5.4) becomes 

— - = F\-x(E hl -E b 1 ) + F x _ 2 (E bl -E h2 ) + F x _ 3 (E bi -E b3 ) ( 5.5 ) 

£ x £ x 
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Table 7 shows the comparison of the hand calculations 
and CFD-ACE software results for wall temperature of 400 F 
and gas temperature of 3500 K between gray surfaces. 


The relative error is due to the fact that in hand 
calculations the outlet is assumed to be isothermal, which 
is not the case for real cases. 



Hand 

Calculation 

Numerical 

Results 

Relative 

Error (%) 

Radiation 

-2512.24 

-2369 

5.701 

to the Wall 




(watts/rad) 





Table 7. Relative Error in Radiative Exchange 


Between Gray Surfaces 


C. GRAY GAS CASE 

A literature search shows that, although there is an 
extensive search for the gray gas approximation, the 
pressure and temperature values are limited to the 
temperatures and pressure values lower than that of 
interest in this study. [2,6,7,8] . Next, CFD-ACE results are 
compared with hand calculations for different optical 
thickness values for parallel plate and round tube cases 
with participating medium. Parallel plate and round tube 
cases are considered consecutively. 


35 




1 . 


Parallel Plate 


In parallel plate case, a gray, nonscattering medium 
is contained between two parallel gray plates. The medium 
is isothermal at temperature T g , with constant absorption 

coefficient k . The two plates are isothermal at temperature 
T w , have the same gray-diffuse emissivity , and are spaced 
a distance D apart as shown in Figure 11. 



Figure 11. Coordinate System and Dimensions for 

the Parallel Plate 

Heat flux at for the specified conditions is given in 
[1] as 


q(T) = cr(T n 4 


T 4 \ 2 [E 3 (t)-E 3 (t l -t)\ 

l + (--l)[l-2 E 3 (t l )\ 
£ 


(5.6) 


At the wall x is equal to zero. So, the equation 
becomes 


q(T = 0) = a(T w - T ;) 


4 m 4 X 2[£ , 3 (0) E 3 (t l )\ 


l + (--l)[l-2 E 3 (t l )\ 
£ 


(5.7) 
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and 


q(r = 0) = g(T j - T g 4 )- [1 E?ATl )] - (5.8) 

l + (—l)[l-2 E 3 (t l )] 

£ 

Equation 5.8 is valid for all the optical thickness 
values other than zero and infinity. For optically thick 
and optically thin cases, two different solutions are 
necessary. 

In optically thick case x L goes to infinity and 
Equation (5.8) becomes 

q = £G(T H 4 -T g 4 ) (5.9) 

In the optically thin situation the following result 
for total heat loss is used for any isothermal volume 
without self absorption in [1] 

Q = 4m 2 oT g 4 V (5.10) 


Multiplying and dividing this equation with 2R and 
writing the heat flux gives 

o 4T,n 2 oT o 4 V/2R 2 T,n 2 oT 4 

„ _ _ _ _ _ 2 _ _ _ _ O / ^ 11 ) 

q - A s ~ (2 m + 2 m + 2HD) ~ 1 D_ D 

H i 


where, V is total volume and A s is the surface area. 

Assumption of £ is greater than D and H is a very good 
one and equation (5.11) becomes 


q = 


4T L n 2 cfiyi2R 2 T L n 2 oT g 


Q__ ___ 

A s (2 Hi + 2DI + 2 HD) } D_ 

H 


(5.12) 
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Since the case is optically thick, taking two values 
for t l is good enough, such as zero and one. The other 
unknown is D/H ratio and for this values, zero and one and 
a half is chosen. Zero corresponds to a very long parallel 
plate, which makes H equal to infinity, but it causes to 
deteriorate the true geometry. On the other hand, one and a 
half corresponds to the real geometry. Figure 12 and 13 
compare the theoretical radiation to the wall values with 
CFD-ACE results for various optical thicknesses, for 
optically thick and optically thin case in a parallel 
plate. 



Figure 12. Comparison of Theoretical and Numerical 

Results for Parallel Plate 
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Figure 13. 


Comparison of Theoretical and Numerical 
Results for Parallel Plate (magnified) 


2. Round Tube 


In round tube case, a gray, nonscattering medium is 
contained between the walls of the cylinder. The medium is 
isothermal at temperature T g , with constant absorption 

coefficient k. The wall of the cylinder is isothermal at 
temperature T w , have the gray-diffuse emissivity , and the 
radius is R as shown in Figure 14. 

The nondimensional heat loss from a gray, 
nonscattering, isothermal cylinder is given in [1] as 


w = 


q(? R ) 


n oT g -J, 


(5.13) 
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Since the walls are not black, but are gray, diffuse 
emitters and reflectors, the heat flux to the wall is no 


longer 7il b , but must be replaced by the radiosity. 


q (T „) = r(n 2 &r; - jj = (j, 

1 - £ 


(5.14) 


Since the walls are not black, but are gray, diffuse 
emitters and reflectors, the heat flux to the wall is no 
longer 7il b , but must be replaced by the radiosity. 


q(T R ) = Win 1 dr; -JJ = (J w ~E b J £ 


l-£„ 


So, the radiosity becomes 


J.. 


(Wn 2 oT 4 ) + ( oT w 4 ) 
1 - 


1 -£ 


+ ¥ 


And the heat flux is 


q(Tj = OW 


0 n% 4 )- 


(¥n 2 T 4 ) + (- £ ^T w 4 ) 
1 - £,„ 


1 -£ 


+ ¥ 


(5.14) 


(5.15) 


(5.16) 


Equation 5.16 is valid for all the optical thickness 
values other than zero and infinity. For optically thick 
and optically thin cases, two different solutions are 
necessary. 

In optically thick case v|/ goes to one according to 
Table 8 and Equation (5.16) becomes 
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r ) — d" 


(5.17) 


(»V)+(r £L! '.‘) 
(»%') -— 


In the optically thin situation the following result 
for total heat loss is used for any isothermal volume 
without self absorption in [1] 

Q = 4m 2 oT g 4 V (5.18) 


Multiplying and dividing this equation with R and 
writing the heat flux gives 


q = 


Q 

A 


4 l R H ^ VIR = 2 T R n 2 oT 4 (—^— 
(ttR 2 +7TR- +2tiRH) s R l + 1 

H 


(5.11) 


where, V is total volume of the cylinder and A s is the 
surface area. 


Since the case is optically thick, taking two values 
for t r is good enough, such as zero and one. The other 
unknown is R/H ratio and for these values, zero and one 
quarter is chosen. Zero corresponds to a very long 
cylinder, which makes H equal to infinity, but it causes to 
deteriorate the true geometry. On the other hand, one 
quarter corresponds to the real geometry. Figure 15 and 16 
compares the theoretical radiation to the wall values with 
CFD-ACE results for various optical thicknesses, for 
optically thick and optically thin case in a round tube. 
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radiation to the wall (W/rad) 



Figure 15. Comparison of Theoretical and Numerical 

Results for Round Tube 



optical thickness 

Figure 16. Comparison of Theoretical and Numerical 

Results for Round Tube (magnified) 


43 















































It is obvious that, as the participating gas gets 
optically thick the results do not agree since, in reality, 
the gas is not isothermal across the chamber cross-section. 

When the gas is optically thick, its absorption 
coefficient is large, so the typical absorption length 
scale is small. In other words, emission from the gas is 
reabsorbed over a very short distance. In such a scenario, 
the radiation on the wall comes from the region of gas 
close to the wall where the gas is at a lower temperature 
than the core temperature (due to the boundary layer nature 
of the temperature profile), and is hence a lower heat flux 
value. 

An important point is to investigate if these hand 
calculations contain radiation to the wall according to the 
expressions used herein or not. So, in order to compare 
these results with each other, the value of radiation to 
the wall for gray walls, transparent medium from CFD-ACE, 
is not added to the hand calculation results. This is shown 
in Appendix A. 

D. NON-GRAY GAS CASE 

A literature search shows that only a limited number 
of non-gray cases have been treated in the past, and have 
been restricted to pure gases at temperatures and pressures 
lower than that of interest in this study [4,5] . The 

radiation heat flux in CO for 1000 °R gas temperature at 1 
atmosphere pressure and 500 °R wall temperature with a 
parabolic initial velocity distribution is considered. The 
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parabolic velocity profile at the inlet is created by using 
the following expression and shown if Figure 17. 


U = 2U m <X- 



(5.12) 



Figure 

Where, 


U 


17 . 

U ave iS 

Re* ,4/ 

ave vfc p. 

p*D 


Parabolic Velocity Profile at the Inlet 
found from the following expression. 

(5.13) 


By using the fire prediction program RADCAL, the 
intensities and transmissivities are found based on the 
wavelength as shown below in Figure 18. 

The numerical result and published result are compared 
in Figure 19 and compared with [4] . The published results 
and numerical results agree well. 
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Rq/kT transmissivity 



8. RADCAL Output for Transmissivity of 

Pure CO 
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Figure 19. 


CFD Output for the Specified Conditions 
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VI. RESULTS AND DISCUSSIONS 


The model problem being treated in this study contains 
several parameters, such as the gas optical thickness, the 
conduction to radiation parameter, the emissivity and 
temperature of the wall, the gas mixture constituents and 
its temperature, etc. Of these, only the parameters of 
primary importance were varied in a parametric study as 
described below so as to determine the role of the main 
contributing factors to the radiative heat loading of the 
duct wall. 

The results presented here are obtained using a 
100x100 grid. Further subdivision of the enclosure does not 
yield significantly more accurate results under the 
specified convergence criteria. Convergence is checked 
using a tolerance of le-10. The CPU time required to 
generate the results varied from 22 minutes to 107 minutes 
(per run) depending on the complexity of the case. 

A. GRAY GAS CASE 

Under the gray gas approximation, a parametric study 
was conducted by varying gas temperature, wall temperature, 
wall emissivity, and optical thickness. The role of the gas 
temperature for the gray gas approximation is shown in 
Figure 20 and 21. Figure 20 shows the convective and 
radiative heat fluxes to the wall versus the optical 
thickness for a low gas to wall temperature ratio and for a 
high gas to wall temperature ratio. It is evident that 
there is an optimal optical thickness at which the 
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radiative flux is a maximum, which can be seen clearly in 
Figure 20. As expected the radiative flux is higher on the 
black wall, and higher gas temperatures. In both figures, 
it is shown that, the radiative flux is a significant 
fraction of the total heat flux especially in the large 
gas-to-wall temp ratio case. 
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Figure 2 0. 


Convective and Radiative Heat Flux for 


1500K Gas 800K Wall Temperature 


The dependence on wall emissivity for the gray gas 
approximation is shown in Figure 22 and Figure 23. These 
figures show the relative contribution of the radiative 
heat flux versus the non-dimensional length along the duct 
for low gas to wall temperature and for high gas to wall 
temperature ratios, respectively. An optical thickness of 3 
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is used which is close to the optimal value as shown in the 
previous figures for maximum radiative flux. It is clear 
that the relative radiative contribution for the high gas 
temperature is almost twice that for the low temperature 
case. 



0 2 4 6 8 10 

optical thickness 

Figure 21. Convective and Radiative Heat Flux for 

1500K Gas 800K Wall Temperature 

Figures 24 and Figure 25 show the radiative heat flux 
versus the non-dimensional length along the duct for 
different wall temperatures at 3500 K gas temperature. 

It is clear in Figure 25 that, the absolute heat flux 
values are nearly the same, although in Figure 24, the 
relative contribution is slightly higher for the hotter 
wall. 
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Figure 22. Relative Contribution of the Radiative 

Heat Flux in Different Wall Emissivities 
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Figure 23. Relative Contribution of the Radiative 

Heat Flux in Different Wall Emissivities 
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Figure 24. Relative Contribution of the Radiative 

Heat Flux in Different Wall Emissivities 
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Figure 25. Relative Contribution of the Radiative 

Heat Flux in Different Wall Emissivities 
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B. 


NONGRAY GAS CASE 


A parametric study was conducted by varying gas 
pressure in addition to the parameters which used in the 
gray gas approximation. Figure 26 shows the typical 
constituents in the gas mixture resulting from the complete 
combustion of RP1 fuel with oxygen as oxidant and an 0/F 
ratio of 3.0. Figure 26 is for the particular gas of a 1500 
K gas mixture at a pressure of 25 atmospheres. 



Figure 26. Typical Mole Fraction of A Gas Mixture 


These 

mole fractions are 

used as 

input 

in RADCAL 

and 

obtain 

spectral radiation 

data as 

shown 

in Figure 

27 


where transmissivity is plotted versus wavelength. The 
bandwidth of each constituent is divided into wavenumbers. 

The Figures 28 and 29 discuss inlet and outlet 
emissivity and show the radiative heat flux to the wall for 
low gas to wall temperature ratio, and for high gas to wall 
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temperature ratio at a chamber pressure of 35 atmosphere. 
The blue line is for black inlet & outlet surfaces, for 
which the radiative heat flux is a maximum. 



Figure 27. Typical Spectral Transmissivity 
Variation of the Mixture 

The purple line is for the case of zero inlet/outlet 
emissivity and shows only the core gas contribution. The 
0.5 emissivity curve is an intermediate result between the 
previous cases. 

Figures 30 and 31 discuss the role of wall emissivity 
by comparing black and gray walls and show the relative 
contribution of the radiative heat flux versus the non- 
dimensional length along the duct for low gas to wall 
temperature ratio and high gas to wall temperature ratio at 
a chamber pressure of 35 atmosphere pressures. 


53 















































It is clear that, the heat flux is about 50% higher to 
black walls, which appears to be fairly independent of the 
actual gas temperature. 



(x/R)/RePr 


Figure 28. Heat Flux for High Gas-to-wall 

Temperature Ratio 
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(x/R)/RePr 


Figure 29. Heat Flux for Low Gas-to-wall 

Temperature Ratio 
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(x/R)/Re Pr 

Figure 30. Heat Flux for High Gas-to-wall 

Temperature Ratio for Black Wall 



1.E-06 1.E-05 1.E-04 
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Figure 31. Heat Flux for High Gas-to-wall 

Temperature Ratio for Gray Wall (emissivity=0.5) 

Figures 32 and 33 discuss role of gas pressure at 3500 
K gas temperature and 400 K wall temperature for the case 
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of a black wall and show the variation of heat flux versus 


non-dimensional length for various chamber pressures. Both 
the absolute and relative values of the radiative heat flux 
do not change appreciably due to the fact that the 
absorption coefficients in the most active bands are 
already around their optimal value, and so much larger 
pressures are needed to alter the radiative heat flux 
significantly. 



1.E-06 1.E-05 1.E-04 

(x/R)/RePr 

Figure 32. Relative Contribution of Heat Flux in 

Black Wall for Various Gas Pressures 



1.E-06 1.E-05 1.E-04 

(x/R)/RePr 

Figure 33. Radiative Heat Flux in Black wall for 

Various Gas Pressures 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A numerical model has been developed for finding the 
radiative heat load on the walls of a combustion chamber. 
Such a model is useful for improving and augmenting the 
predictive capabilities of other numerical codes that 
provide only convective heat loads. 

Although some of the trends in results are expected, 
this study provides concrete quantitative estimates of the 
role of gas temperature, wall temperature and emissivity, 
etc. 

It was found that the radiative fluxes are least 
sensitive to the wall temperature and chamber pressure in 
the range considered. However, the effective emissivities 
of the inlet and outlet surfaces were found to be crucial 
in determining the magnitude of these fluxes. 

Most importantly, it was found that, a proper non-gray 
spectral treatment of the gas mixture is critical for an 
accurate analysis. A gray treatment is too simplistic and 
yields unreliable results. 

More accurate and realistic estimates can be obtained 
by expanding the study to include the following: 

• The absorption-emission and scattering effects of 
particulate matter such as soot 

• The role of the liquid fuel that is injected to form a 
film on the chamber walls 
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• The interaction of the combustion mechanics with the 


flow and heat transfer characteristics 

• The effect of the resulting non-uniform distribution 
of gaseous species in the chamber 
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APPENDIX A. THE WALL CONTRIBUTION TO HEAT FLUX 


If the medium is not at radiative equilibrium (i.e. 
conduction and/or convection is not negligible and 
radiative heat transfer is not dominant) the heat flux is 
found by using the definition of the exponential integrals 
in [1] as 

t n 

q(T) = 2 J ] E 3 (t) - 2 J 2 E 3 (t l -T) + 2 oT iu 4 \E 2 (t- T)dr'- 2 oT m 4 J E 2 (t -r)dr' 

0 T 

and it becomes 

q(T) = 2J l E 3 (T)-2J 2 E 3 (r L -T) + 2oT m 4 [E 3 (T)-E 3 (T L - T)] 

l 



Figure 34. Optical Coordinates for the Model 

As shown in Figure 34 optical coordinate at the upper 
wall is zero and optical coordinate at the lower wall is x L . 

Setting x equal to zero and x L in the previous equation 
consecutively gives 

<7(0) = 2J x E 3 (0) - 2 J 2 E 3 (t l ) + 2oT m 4 [E 3 (0) - E 3 (t l )] = (cfT, 4 - J ,) 

l-^i 

q(r L ) = 2J\E 3 (t l ) -2 J 2 E 3 (0) + 2ctT m 4 [E 3 (r) - E 3 (0)] = (oT 2 4 - J 2 ) 

\-e 2 
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As seen in the equation above, the expression found 
for the heat flux at the upper wall considers the 
emissivity and temperatures of the top and bottom walls. 
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APPENDIX B. TYPICAL CFD-ACE OUTPUT 
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: Windows 
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Zone No. 


1 


VC Name. : NoName 

No. of Cells : 9801 

Material Type : Fluid 

Density method : Constant = 2.97E+00 

Viscosity method : Constant_Dyn = 9.52E-05 

Cond. method: Prandtl = 6.5E-01 

Sp. Heat method: Constant = 1.94E+03 

Absorption Coeff. Set: gas 
Emissivity Set: None 


RAD_SNORD Statistics 


Parameters 


Radiation 
Gray model 
CASN-DOM 

S4 scheme 
sn_hrange scheme 
% Relative Error 
err_0m 
err_lm_mu 
err_lm_xi 
err_lm_et 
err_2m_mu 
err_2m_xi 
err_2m_et 


= TRUE 

= TRUE 

= TRUE 

= TRUE 
in Moments 
-2.980232E-06 
-1.092752E-05 
-1.092752E-05 
-1.000000E+02 
-1.185405E-05 
-1.185405E-05 
-1.185405E-05 
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Total Mass Flow Summary 


2.786166E-02 


2.786166E-02 -5.55112E-17 
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Name 


Key Type 


X-axis 


Y-axis 


Z-axis 

1 NoName 28 Wall 0.000000E+00 0.000000E+00 

1.778739E+02 


Viscous Moments 


Name 


Z-axis 


1 NoName 
-1.392094E-03 


Key Type X-axis Y-axis 


28 Wall 0.000000E+00 0.000000E+00 


Boundary-by-Boundary Heat Transfer Summary (watts/rad) 


Name 

key 

Type 

COND.+CONV. 

W_SRC.+CVD 

RADIATION* 

Sum 

NoName 

37 

Inlet 

1.8899E+05 

0.0000E+00 

1.2769E+02 


1.8912E+05 

NoName 

36 

Outlet 

-1.7 693E+05 

0.0000E+00 

3.6525E+02 - 


1.7656E+05 

NoName 

28 

Wall 

-8.2173E+03 

0.0000E+00 

-4.3441E+03 - 



1.2561E+04 


Total volume source 
0.0000E+00 
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Total Heat Imbalance 3.8479E+03 0.0000E+00 -3.8511E+03 - 

3.2126E+00 


Total wall HEAT source 
0.0000E+00 

* Radiation Summary Convention is +ve if emmitting, -ve if absorbing. 

* Cond+Conv Summary Convention is +ve if flux into the cell, -ve if 
flux leavi ng the cell. 


End of Iterative Cycle 


Final Time Elapsed Time= 6.552923E+02 Delta-time= 6.552923E+02 
Normal Termination 
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APPENDIX C. TYPICAL RADCAL INPUT 


l 


0.2850 

1000. 7.5 15.6 

. 0 

. 0 

. 0 

77.0 

. 0 

0. 50. 

1 

10000 . 






1.000 

1000. 7.5 15.6 . 

.0 , 

.0 , 

. 0 

77.0 , 

. 0 

0. 2000. 3000. 






0 

0.030 

1000. 7.5 15.6 . 

.0 , 

.0 , 

. 0 

77.0 , 

. 0 

1500 . 

50. 10000. 






1 

0.010 

1000. 7.5 15.6 . 

.0 , 

.0 , 

. 0 

77.0 . 

. 0 

1500 . 

50. 10000. 






1 

0.003 

1000. 7.5 15.6 . 

.0 , 

.0 , 

. 0 

77.0 , 

. 0 

1500 . 

50. 10000. 






1 

0.001 

1000. 7.5 15.6 . 

.0 , 

.0 , 

. 0 

77.0 , 

. 0 

1500 . 

50. 10000. 






1 

0.0003 

1000. 7.5 15.6 

. 0 

. 0 

. 0 

77.0 

. 0 

1500 . 

50. 10000. 






1 

0.0001 

1000. 7.5 15.6 

. 0 

. 0 

. 0 

77.0 

. 0 


1500. 50. 10000. 
0 


Data input file for RADCAL 

line 1: number of homogeneous elements, n 

line 2: pathlength (m), temperature (K), C02 (kPa), H20, CH4, CO, 

02, N2, fv 

lines 3 through n+1: same as line 2 for the rest of the elements 
line n+2: wall temperature (K), minimum wavenumber (cm-1), maximum 
wavenumber 

line n+3: 0, or the number of homogeneous elements in the next case 
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APPENDIX D. TYPICAL RADCAL OUTPUT 


Radial Profiles 


Partial Pressures, kPa 

J dist,m temp,K C02 H20 CH4 CO 

1 .2850 1000. 7.500 15.600 .000 .000 

wall 400. 


02 N2 FV 
.000 77.0 0.0 


Total directional radiated energy flux = .307139E+04 Watts/m-2/strad 


Spectral Intensity Distribution, Watts/m-2/micron/strad 


micron 

intensity 

tau 

micron 

intensity 

tau 

1.005 

.0000E+00 

1.0000 

6.061 

.2463E+03 

. 8353 

1.010 

. 0000E+00 

1.0000 

6.154 

.1477E+03 

.8975 

1.015 

. 0000E+00 

1.0000 

6.250 

.1126E+03 

. 9189 

1.020 

.0000E+00 

1.0000 

6.349 

. 1516E+03 

. 8865 

1.026 

. 0000E+00 

1.0000 

6.452 

.3674E+03 

. 7137 

1.031 

. 0000E+00 

1.0000 

6.557 

.5144E+03 

.5824 

1.036 

. 0000E+00 

1.0000 

6.667 

. 4936E+03 

.5821 

1.042 

.0000E+00 

1.0000 

6.780 

. 4175E+03 

. 6309 

1.047 

. 0000E+00 

1.0000 

6.897 

. 3891E+03 

. 6404 

1.053 

.0000E+00 

1.0000 

7.018 

.3447E+03 

. 6665 

1.058 

. 0000E+00 

1.0000 

7.143 

.3104E+03 

. 6852 

1.064 

.0000E+00 

1.0000 

7.273 

.2909E+03 

. 6904 

1.070 

. 0000E+00 

1.0000 

7.407 

.2773E+03 

. 6897 

1.075 

. 0000E+00 

1.0000 

7.547 

.2674E+03 

. 6850 

1.081 

.2750E-03 

1.0000 

7.692 

. 2309E+03 

. 7133 

1.087 

. 5450E-03 

1.0000 

7.843 

.1981E+03 

. 7403 

1.093 

. 5559E-02 

1.0000 

8.000 

.1503E+03 

.7915 

1.099 

. 7908E-02 

. 9999 

8.163 

. 1068E+03 

. 8431 

1.105 

. 1017E-01 

.9999 

8.333 

. 7042E+02 

.8902 

1.111 

.1929E-01 

. 9999 

8.511 

. 4206E+02 

. 9303 

1.117 

.3087E-01 

. 9998 

8.696 

.2218E+02 

. 9608 

1.124 

.3884E-01 

. 9998 

8.889 

.1643E+02 

. 9690 

1. 130 

.4722E-01 

. 9998 

9.091 

.2560E+02 

. 9484 

1.136 

.5675E-01 

. 9997 

9.132 

.2681E+02 

. 9452 

1. 143 

.6802E-01 

. 9997 

9.174 

.2729E+02 

. 9434 

1.149 

.7338E-01 

. 9997 

9.217 

.2677E+02 

. 9437 

1. 156 

.9351E-01 

. 9996 

9.259 

.2505E+02 

.9466 

1.163 

.1019E+00 

. 9996 

9.302 

.2200E+02 

. 9524 

1. 170 

. 1129E+00 

. 9995 

9.346 

.17 68E+02 

. 9612 

1.176 

.1128E+00 

. 9996 

9.390 

.1214E+02 

. 9730 

1. 183 

.1404E+00 

. 9995 

9.434 

.1413E+02 

. 9681 

1.190 

.1451E+00 

. 9995 

9.479 

.1813E+02 

. 9585 

1.198 

.9713E-01 

. 9997 

9.524 

.2072E+02 

. 9519 

1.205 

.1105E+00 

. 9996 

9.569 

.2199E+02 

. 9482 
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1.212 

.9183E-01 

. 9997 

9.615 

.2207E+02 

. 9472 

1.220 

.1101E+00 

. 9997 

9.662 

.2119E+02 

. 9486 

1.227 

.5030E-01 

.9999 

9.709 

.1966E+02 

. 9515 

1.235 

.4276E-01 

. 9999 

9.756 

.1780E+02 

. 9555 

1.242 

.4418E-01 

. 9999 

9.804 

.1598E+02 

. 9594 

1.250 

. 1951E-01 

1.0000 

9.852 

.1440E+02 

. 9629 

1.258 

.1816E-01 

1.0000 

9.901 

.1324E+02 

. 9653 

1.266 

.2371E-01 

. 9999 

9.950 

.1255E+02 

.9666 

1.274 

.3800E-01 

. 9999 

10.000 

. 1225E+02 

. 9669 

1.282 

. 5689E-01 

. 9999 

10.050 

.1224E+02 

.9664 

1.290 

.7012E-01 

. 9999 

10.101 

.1226E+02 

. 9658 

1.299 

. 3162E+00 

. 9994 

10.152 

.1212E+02 

. 9657 

1.307 

.3513E+00 

. 9993 

10.204 

.1167E+02 

.9664 

1.316 

.8 4 66E+00 

. 9984 

10.256 

.1080E+02 

. 9684 

1.325 

.2019E+01 

. 9964 

10.309 

.9539E+01 

. 9717 

1.333 

. 2680E+01 

. 9954 

10.363 

.7912E+01 

. 9761 

1.342 

.7485E+01 

. 9876 

10.417 

.6633E+01 

. 9796 

1.351 

. 1995E+02 

. 9682 

10.471 

.8154E+01 

. 9746 

1.361 

.3008E+02 

. 9539 

10.526 

.9307E+01 

. 9705 

1.370 

.3041E+02 

. 9551 

10.582 

.1006E+02 

.967 6 

1.379 

.2758E+02 

. 9608 

10.638 

.1035E+02 

. 9660 

1.389 

.3325E+02 

. 9545 

10.695 

.1024E+02 

. 9658 

1.399 

.3583E+02 

. 9527 

10.753 

.9810E+01 

.9667 

1.408 

.3438E+02 

. 9563 

10.811 

.9186E+01 

. 9683 

1.418 

.3388E+02 

. 9585 

10.870 

.8521E+01 

. 9701 

1.429 

.30 91E+02 

. 9635 

10.929 

. 7865E+01 

. 9719 

1.439 

.2601E+02 

. 9703 

10.989 

.7284E+01 

. 9735 

1.449 

.2227E+02 

. 9755 

11.050 

. 6814E+01 

. 9748 

1.460 

. 1790E+02 

. 9810 

11. Ill 

. 6460E+01 

. 9756 

1.471 

. 1061E+02 

. 9891 

11.173 

. 6270E+01 

. 9759 

1.481 

.9673E+01 

. 9904 

11.236 

.6159E+01 

. 9759 

1.493 

.8050E+01 

. 9923 

11.299 

.6104E+01 

. 9757 

1.504 

. 6925E+01 

. 9936 

11.364 

. 6029E+01 

. 9756 

1.515 

.5957E+01 

. 9947 

11.429 

.6057E+01 

. 9750 

1.527 

.4139E+01 

. 9964 

11.494 

.6112E+01 

. 9743 

1.538 

.3552E+01 

. 9970 

11.561 

. 6162E+01 

. 9736 

1.550 

.2128E+01 

. 9983 

11.628 

. 6207E+01 

. 9729 

1.563 

.1522E+01 

. 9988 

11.696 

. 6248E+01 

. 9722 

1.575 

. 1290E+01 

. 9990 

11.765 

. 6283E+01 

. 9715 

1.587 

.1014E+01 

. 9993 

11.834 

. 6376E+01 

. 9705 

1.600 

.6685E+00 

. 9995 

11.905 

. 6463E+01 

. 9695 

1.613 

.4852E+00 

. 9997 

11.976 

.6542E+01 

. 9685 

1.626 

. 4411E+00 

. 9997 

12.048 

.7069E+01 

. 9653 

1.639 

.4084E+00 

. 9997 

12.121 

.7273E+01 

. 9636 

1.653 

. 3866E+00 

. 9998 

12.195 

. 9147E+01 

. 9533 

1.667 

.3688E+00 

. 9998 

12.270 

. 8254E+01 

. 9570 

1.681 

. 4618E+00 

. 9997 

12.346 

.8842E+01 

. 9530 

1.695 

.5073E+00 

. 9997 

12.422 

. 9667E+01 

. 9475 

1.709 

.6027E+00 

. 9997 

12.500 

. 1025E+02 

. 9432 

1.724 

. 1392E+01 

. 9993 

12.579 

.1104E+02 

. 9376 

1.739 

. 3258E+01 

. 9983 

12.658 

.1197E+02 

. 9309 

1.754 

.9787E+01 

. 9950 

12.739 

.1348E+02 

. 9205 

1.770 

. 2302E+02 

. 9886 

12.821 

.1497E+02 

. 9098 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


786 

.4804E+02 

. 9769 

12.903 

.1582E+02 

. 9026 

802 

.8527E+02 

. 9601 

12.987 

.17 67E+02 

. 8888 

818 

.1165E+03 

. 9469 

13.072 

.1978E+02 

. 8728 

835 

. 1231E+03 

. 9453 

13.158 

.2243E+02 

. 8526 

852 

.1330E+03 

. 9424 

13.245 

.2550E+02 

. 8286 

869 

. 1840E+03 

. 9224 

13.333 

. 2466E+02 

. 8305 

887 

.1082E+03 

. 9555 

13.423 

.3142E+02 

. 7792 

905 

.8995E+02 

. 9639 

13.514 

.38 64E+02 

. 7222 

923 

.1072E+03 

. 9580 

13.605 

.3805E+02 

. 7201 

942 

.1342E+03 

. 9487 

13.699 

.3886E+02 

. 7074 

961 

.1511E+03 

. 9435 

13.793 

. 4769E+02 

. 6325 

980 

.1217E+03 

. 9555 

13.889 

.7088E+02 

.4409 

000 

.1395E+03 

. 9501 

13.986 

.6696E+02 

. 4591 

010 

.9192E+02 

. 9675 

14.085 

.6327E+02 

.4766 

020 

.8857E+02 

. 9690 

14.184 

.6714E+02 

.4310 

030 

.4371E+02 

. 9849 

14.286 

.6781E+02 

.4113 

041 

. 4465E+02 

. 9847 

14.388 

.7217E+02 

.3579 

051 

.4845E+02 

. 9836 

14.493 

. 7992E+02 

.2711 

062 

.3829E+02 

. 9872 

14.599 

.8159E+02 

.2371 

073 

.4135E+02 

. 9863 

14.706 

.7775E+02 

.2545 

083 

.1168E+02 

. 9962 

14.815 

. 8612E+02 

. 1531 

094 

. 4225E+01 

. 9986 

14.925 

.9416E+02 

. 0500 

105 

.1757E+01 

. 9994 

15.038 

.90 91E+02 

. 0589 

116 

.1581E+01 

. 9995 

15.152 

.6921E+02 

.2646 

128 

.1470E+01 

. 9995 

15.267 

.6953E+02 

.2416 

139 

.1308E+01 

. 9996 

15.385 

.6904E+02 

.2267 

151 

. 1150E+01 

. 9996 

15.504 

.5944E+02 

.3162 

162 

.1113E+01 

. 9997 

15.625 

.5130E+02 

.3937 

174 

.1056E+01 

. 9997 

15.748 

.4844E+02 

.4116 

186 

. 1081E+01 

. 9997 

15.873 

.4464E+02 

.4427 

198 

. 1131E+01 

. 9997 

16.000 

.3964E+02 

.4913 

210 

. 1246E+01 

. 9996 

16.129 

.3698E+02 

.5119 

222 

.1342E+01 

. 9996 

16.260 

.4188E+02 

. 4313 

235 

.1390E+01 

. 9996 

16.393 

.3814E+02 

.4672 

247 

. 2010E+01 

. 9994 

16.529 

.3307E+02 

.5245 

260 

.2084E+01 

. 9994 

16.667 

.3046E+02 

.5491 

273 

. 2685E+01 

. 9992 

16.807 

.3175E+02 

.5159 

286 

.3155E+01 

. 9991 

16.949 

.2961E+02 

.5350 

299 

.4134E+01 

. 9988 

17.094 

.2789E+02 

.5486 

312 

.4731E+01 

. 9987 

17.241 

. 2691E+02 

.5512 

326 

.6020E+01 

. 9983 

17.391 

. 2615E+02 

.5503 

339 

.6801E+01 

. 9981 

17.544 

. 2582E+02 

. 5421 

353 

. 8293E+01 

. 9977 

17.699 

. 2576E+02 

.5287 

367 

. 9228E+01 

. 9975 

17.857 

. 2516E+02 

.5249 

381 

.1570E+02 

. 9958 

18.018 

. 2490E+02 

. 5147 

395 

.1983E+02 

. 9947 

18.182 

.2457E+02 

.5054 

410 

.3383E+02 

. 9910 

18.349 

.2429E+02 

.4950 

424 

.52 94E+02 

. 9860 

18.519 

. 2388E+02 

.4869 

439 

.9285E+02 

. 9755 

18.692 

. 2350E+02 

.4780 

454 

.1579E+03 

. 9586 

18.868 

. 22 98E+02 

.4722 

469 

. 2717E+03 

. 9291 

19.048 

. 2256E+02 

.4641 

484 

.4305E+03 

. 8883 

19.231 

.2215E+02 

. 4556 

500 

.6355E+03 

. 8359 

19.417 

.2176E+02 

.4462 
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2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


516 

. 8206E+03 

.7891 

19.608 

.2134E+02 

. 4374 

532 

. 9776E+03 

.7500 

19.802 

.2090E+02 

. 4293 

548 

. 1067E+04 

.7282 

20.000 

.2043E+02 

. 4217 

564 

.9858E+03 

.7500 

20.202 

.1998E+02 

.4136 

581 

.1036E+04 

.7384 

20.408 

. 1952E+02 

.4060 

597 

. 7898E+03 

. 8012 

20.619 

. 1903E+02 

.3989 

614 

. 7146E+03 

. 8208 

20.833 

.1853E+02 

.3923 

632 

.5482E+03 

.8630 

21.053 

.1803E+02 

.3861 

64 9 

.4264E+03 

.8938 

21.277 

.1755E+02 

.3790 

667 

.1657E+04 

. 5883 

21.505 

.1706E+02 

. 3724 

685 

.1849E+04 

.5419 

21.739 

. 1657E+02 

.3662 

703 

.1975E+04 

. 5117 

21.978 

.1607E+02 

.3603 

721 

.1645E+04 

.5943 

22.222 

.1557E+02 

. 3549 

740 

.1340E+04 

. 6702 

22.472 

. 1506E+02 

.3499 

759 

. 1206E+04 

.7038 

22.727 

.1456E+02 

.3452 

778 

. 1376E+04 

. 6625 

22.989 

.1406E+02 

.3408 

797 

.1440E+04 

. 6473 

23.256 

.1356E+02 

.3366 

817 

.1471E+04 

. 6400 

23.529 

.1307E+02 

.3326 

837 

.1373E+04 

.6642 

23.810 

. 1258E+02 

. 3291 

857 

.1190E+04 

.7092 

24.096 

. 1210E+02 

.3258 

878 

. 9948E+03 

.7570 

24.390 

.1163E+02 

.3226 

899 

.8855E+03 

.7837 

24.691 

.1116E+02 

. 3196 

920 

. 7522E+03 

.8163 

25.000 

.1071E+02 

.3168 

941 

.64 4 6E+0 3 

. 8425 

25.316 

.1025E+02 

.3150 

963 

. 5486E+03 

.8659 

25.641 

.9802E+01 

.3134 

985 

. 4688E+03 

. 8853 

25.974 

.9368E+01 

.3117 

008 

.3842E+03 

. 9059 

26.316 

. 8945E+01 

.3102 

030 

.3267E+03 

. 9198 

2 6.667 

.8535E+01 

.3088 

053 

. 2719E+03 

. 9332 

27.027 

. 8122E+01 

.3086 

077 

. 2190E+03 

. 9461 

27.397 

.7724E+01 

.3085 

101 

.1660E+03 

. 9590 

27.778 

. 7339E+01 

. 3084 

125 

. 1280E+03 

. 9683 

28.169 

.6967E+01 

.3083 

150 

.1127E+03 

. 9720 

28.571 

.6610E+01 

. 3082 

175 

.9436E+02 

. 9765 

28.986 

.6254E+01 

.3094 

200 

.7517E+02 

. 9812 

29.412 

.5913E+01 

. 3105 

226 

.6378E+02 

. 9840 

29.851 

.5585E+01 

.3116 

252 

.5180E+02 

. 9870 

30.303 

.5271E+01 

.3127 

279 

.4519E+02 

. 9886 

30.769 

. 4971E+01 

.3137 

306 

.3939E+02 

. 9900 

31.250 

. 4673E+01 

.3162 

333 

.3259E+02 

. 9917 

31.746 

. 4389E+01 

. 3186 

361 

.2895E+02 

. 9926 

32.258 

. 4118E+01 

.3210 

390 

. 2597E+02 

. 9933 

32.787 

. 3860E+01 

.3234 

419 

.2072E+02 

. 9946 

33.333 

.3615E+01 

. 3258 

448 

.17 64E+02 

. 9954 

33.898 

.3374E+01 

.3295 

478 

. 1464E+02 

. 9961 

34.483 

.3145E+01 

.3333 

509 

. 1238E+02 

. 9967 

35.088 

.2928E+01 

. 3371 

540 

. 1063E+02 

. 9972 

35.714 

. 2723E+01 

.3408 

571 

. 8686E+01 

. 9977 

36.364 

. 2529E+01 

.3445 

604 

.6410E+01 

. 9983 

37.037 

. 2338E+01 

.3505 

636 

.4385E+01 

. 9988 

37.736 

.2158E+01 

.3565 

670 

.3877E+01 

. 9989 

38.462 

. 1988E+01 

.3626 

704 

. 3277E+01 

. 9991 

39.216 

. 1829E+01 

.3687 

738 

. 2320E+01 

. 9993 

40.000 

. 1679E+01 

. 3749 
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3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 


774 

. 2023E+01 

. 9994 

40.816 

.1534E+01 

.3832 

810 

. 1604E+01 

. 9995 

41.667 

. 1398E+01 

.3918 

846 

.1533E+01 

. 9996 

42.553 

.1272E+01 

.4005 

883 

. 1524E+01 

. 9996 

43.478 

.1154E+01 

.4095 

922 

.1518E+01 

. 9995 

44.444 

.1044E+01 

.4187 

960 

. 1460E+01 

. 9996 

45.455 

.9412E+00 

. 4290 

000 

.1480E+01 

. 9995 

46.512 

.8458E+00 

.4395 

040 

. 1690E+01 

. 9995 

47.619 

. 7576E+00 

. 4505 

082 

.1985E+01 

. 9994 

48.780 

.67 63E+00 

.4619 

124 

.2304E+01 

. 9993 

50.000 

.6015E+00 

.4737 

167 

. 2646E+01 

. 9991 

51.282 

.5334E+00 

.4855 

211 

. 2953E+04 

. 0324 

52.632 

. 4711E+00 

.4978 

255 

. 2995E+04 

. 0031 

54.054 

.4141E+00 

.5106 

301 

. 2954E+04 

. 0008 

55.556 

. 3623E+00 

. 5241 

348 

. 2905E+04 

. 0006 

57.143 

.3152E+00 

.5383 

396 

. 2846E+04 

. 0039 

58.824 

.2731E+00 

.5525 

444 

. 2625E+04 

.0648 

60.606 

. 2351E+00 

.5675 

494 

. 1643E+04 

. 4038 

62.500 

. 2009E+00 

.5835 

545 

. 57 66E+03 

.7868 

64.516 

.1703E+00 

. 6006 

598 

.1604E+03 

. 9395 

66.667 

.1430E+00 

. 6190 

651 

.7032E+02 

. 9729 

68.966 

.1200E+00 

. 6354 

706 

.6562E+02 

. 9742 

71.429 

.9956E-01 

. 6531 

762 

.8193E+02 

. 9671 

74.074 

.8162E-01 

. 6723 

819 

.9943E+02 

. 9592 

76.923 

.6593E-01 

. 6933 

878 

.1257E+03 

. 9472 

80.000 

.5227E-01 

. 7166 

938 

.1503E+03 

. 9354 

83.333 

.4202E-01 

. 7328 

000 

.1674E+03 

. 9263 

86.957 

.3324E-01 

. 7503 

063 

.1852E+03 

. 9164 

90.909 

. 2578E-01 

.7695 

128 

.1905E+03 

. 9118 

95.238 

. 1949E-01 

.7909 

195 

. 2603E+03 

. 8764 

100.000 

.1424E-01 

. 8150 

263 

.2817E+03 

.8625 

105.263 

.1073E-01 

. 8295 

333 

. 3009E+03 

. 8491 

111.111 

.7875E-02 

. 8452 

405 

.3575E+03 

.8155 

117.647 

.5586E-02 

.8625 

479 

.3563E+03 

.8106 

125.000 

.3781E-02 

. 8818 

556 

. 4000E+03 

.7808 

133.333 

.2387E-02 

. 9038 

634 

.4012E+03 

. 7733 

142.857 

.1644E-02 

. 9130 

714 

. 4043E+03 

.7642 

153.846 

.1088E-02 

. 9228 

797 

. 4231E+03 

.7450 

166.667 

.6855E-03 

. 9333 

882 

.3710E+03 

.7687 

181.818 

.4039E-03 

. 9445 

970 

. 2889E+03 

.8135 

200.000 

. 2162E-03 

. 9567 


The effective absorption coef. is .654587E+00/m 
The Planck-mean absorption coef. is .290665E+01/m 
The wall-incident mean is .290665E+01/m 


Radial Profiles 


Partial Pressures, kPa 
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J dist,m temp,K C02 H20 CH4 CO 02 N2 

FV 

1 1.0000 1000. 7.500 15.600 .000 .000 .000 77.000 

.0000E+00 

wall 0. 

Total directional radiated energy flux = .125184E+04 Watts/m-2/strad 


Spectral Intensity Distribution, Watts/m-2/micron/strad 


micron 

intensity 

tau 

micron 

intensity 

tau 

3.339 

.1085E+03 

. 9723 

4.008 

.5337E+01 

. 9984 

3.367 

.9697E+02 

. 9751 

4.049 

.6130E+01 

. 9981 

3.396 

.8553E+02 

. 9779 

4.090 

.7178E+01 

. 9977 

3.425 

.6933E+02 

. 9820 

4.132 

.8308E+01 

. 9974 

3.454 

.58 94E+02 

. 9846 

4.175 

.2151E+04 

. 3038 

3.484 

.4918E+02 

. 9870 

4.219 

.3041E+04 

. 0006 

3.515 

.4178E+02 

. 9889 

4.264 

.2995E+04 

. 0000 

3.546 

.3564E+02 

. 9905 

4.310 

.2947E+04 

. 0000 

3.578 

.2868E+02 

. 9923 

4.357 

.2897E+04 

. 0000 

3.610 

.2096E+02 

. 9943 

4.405 

.2847E+04 

. 0000 

3.643 

.1498E+02 

. 9959 

4.454 

.27 93E+04 

. 0012 

3.676 

.1314E+02 

. 9964 

4.505 

.2490E+04 

.0932 

3.711 

.1080E+02 

. 9970 

4.556 

.1262E+04 

.5315 

3.745 

.7920E+01 

. 9978 

4.608 

.4181E+03 

. 8417 

3.781 

.67 93E+01 

. 9981 

4.662 

.2144E+03 

. 9172 

3.817 

.5574E+01 

. 9984 

4.717 

. 2115E+03 

. 9166 

3.854 

.5368E+01 

. 9984 

4.773 

. 2523E+03 

.8983 

3.891 

.5338E+01 

. 9984 

4.831 

. 2964E+03 

. 8778 

3.929 

.5280E+01 

. 9984 

4.890 

.3505E+03 

. 8522 

3.968 

.5131E+01 

. 9985 

4.950 

.3933E+03 

. 8302 


The effective absorption coef. is .718888E-01/m 
The Planck-mean absorption coef. is .199213E+01/m 
The wall-incident mean is .199213E+01/m 
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APPENDIX E. TYPICAL TEP OUTPUT 


INPUT FILE:D:\PC_TEP\tep.dat 
OUPUT FILE:D:\PC_TEP\tep.out 

THERMODYNAMIC DATA FILE:D:\PC_TEP\species\chno.dat 
TITLE RPI/G02 
DATA 
$DATA 
ODE=l , 

NZONES=l , 

ASUB(1)=2, 

NASUB=1, 

ASUP(1)= 6, 

NASUP=1, 

$END 

REACTANTS 

298.15F .773 

298.150 

NAMELISTS 

$ODE 

RKT = .TRUE. 

PSIA=F, 

T(1)=1500, 

P (1)=25, 

OFSKED(1)=3, 

OF= T, 

$END 

OTITLE RPI/G02 
ODATA 

$DATA 
ODE=l, 

NZONES=l, 

ASUB(1)=2, 

NASUB=1, 

ASUP(1)= 6, 

NASUP=1, 

$END 


C 1. 
0 2 . 


H 1.9423 


100 . 

100 . 


-5430.L 
0.0G 


************************************************************************************************** 

****************** 

CALCULATE ODE AREA RATIO AND PRESSURE SCHEDULES FOR ZONE 1 

************************************************************************************************** 

****************** 

REACTANTS 


C 1.0000 H 

0.77300 

1.9423 

0.0000 

0.0000 

0.0000 

100.000000 

-5430.00 

L 

298.150 

O 2.0000 

0.0000 

0.0000 

0.0000 

0.0000 

100.000000 

0.00 

G 

298.150 


0.00000 

NAMELISTS 

$ODE 

RKT = .TRUE. 

PSIA=F, 

T(1)=1500, 

P(1)=25, 

OFSKED(1)=3, 

OF= T, 

$END 

0SPECIES BEING CONSIDERED IN THIS SYSTEM 


J 3/78 
9/65 C02 

C 

J12/67 

CH 

L 5/84 

CH4 

J 9/65 

CO 

J12/69 
6/88 JET- 

C2 

-A(G) 

J 3/61 

C2H2 

J12/69 

C3 

J12/69 

C4 

DR 9 

3/77 H2 

C12H26 

J 3/77 

H 

J12/70 

HCO 

J 9/78 

H02 
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J 3/79 H20 

J 3/77 02 

J 3/78 C(GR) 

PI 9 C12H2 6(L) 

OOF = 3.000000 

ENTHALPY 

(KG-MOL)(DEG K)/KG 
OKG-ATOMS/KG 
C 
H 
0 

ENTHALPY IN BTU/LBM : 
FROM REACTANTS : 
FROM DELH( ) : 

FROM DELHI( ) 
TOTAL : 

1 

EXPANSION 


L 6/80 H202 

J 3/79 H20(L) 


EFFECTIVE FUEL 
HPP (2) 

-0.19573987E+03 
BOP (1,2) 

0.71587443E-01 
0.13904428E+00 
0.00000000E+00 


-174.9239 
0.0000 
0.0000 
-174.9239 

THEORETICAL ROCKET 


J 3/77 0 

Est. JP10(L) 


EFFECTIVE OXIDANT 
HPP (1) 

0.00000000E+00 
BOP (1,1) 
0.00000000E+00 
0.00000000E+00 
0.62502339E-01 


J 6/77 OH 
L 6/88 JET-A(L) 


MIXTURE 

HSUBO 

-0.48934967E+02 
BO (I) 

0.17896861E-01 
0.34761071E-01 
0.46876755E-01 


COMPOSITION DURING 


ZONE = 1 

PERFORMANCE ASSUMING EQUILIBRIUM 


OPC = 367.4 PSIA 


FROM AN ASSIGNED PRESSURE AND TEMPERATURE 

WT FRACTION ENTHALPY 


STATE TEMP DENSITY 

CHEMICAL FORMULA (SEE NOTE) CAL/MOL 

DEG K G/CC 

FUEL C 1.00000 H 1.94230 1.00000 -5430.000 

L 298.15 0.7730 

OXIDANT O 2.00000 1.00000 0.000 

G 298.15 0.0000 

0O/F=3.OOOOE+OO PERCENT FUEL=2.5000E+01 EQUIVALENCE RATIO=l.1343E+00 STOIC MIXTURE RATIO=3.4030E+00 
DENSITY=0.0000E+00 


1 

CHAMBER 

THROAT 

EXIT 

EXIT 

PC/P 

1.0000 

1.7704 

1.0598 

39.796 

P, ATM 

25.00 

14.12 

23.59 

0.6282 

T, DEG K 

1500 

1365 

1486 

822 

H, CAL/G 

-1845.1 

-1902.4 

-1851.2 

-2135.1 

S, CAL/ (G) (K) 

2.0645 

2.0645 

2.0645 

2.0645 

G, CAL/GRAM 

-4941.9 

-4721.1 

-4918.6 

-3832.8 

U, CAL/GRAM 

-1950.2 

-1998.1 

-1955.3 

-2192.6 

DEN (G/LITER) 

5.76E+00 

3.57E+00 

5.48E+00 

2.65E-01 

M, MOL WT 

28.347 

28.347 

28.347 

28.427 

(DLV/DLP)T 

- 1.00000 

- 1.00000 

- 1.00000 

-1.00502 

(DLV/DLT)P 

1.0000 

1.0000 

1.0000 

1.0833 

CP, CAL/(G)(K) 

0.4269 

0.4239 

0.4265 

0.5189 

CP GAS(SF) 

0.4137 

0.4057 

0.4129 

0.3650 

GAMMA GAS(SF) 

1.2039 

1.2087 

1.2043 

1.2367 

GAMMA (S) 

1.1963 

1.1980 

1.1965 

1.1806 

SON VEL,M/SEC 

725.5 

692.6 

722.1 

532.9 

MU, POISE 

5.24E-04 

4.89E-04 

5.20E-04 

3.32E-04 

K,ERG/S-CM-K 

1.41E+04 

1.30E+04 

1.40E+04 

8.75E+03 

PRANDTL NO 

0.64509 

0.63850 

0.64448 

0.57932 

MACH NUMBER 

0.0000 

1.0000 

0.3124 

2.9246 

AE/AT 


1.0000 

2.0000 

6.0001 

CSTAR, FT/SEC 


3358 

3358 

3358 

CF VAC 


1.242 


1.673 

CF 


0.677 


1.522 

IVAC,LBF-S/LBM 


129.58 


174.65 

I, LBF-SEC/LBM 


70.63 


158.91 

MOL WT(MIX) 

28.347 

28.347 

28.347 

28.427 

MOLE FRACTIONS 





CH4 

0.000000 

0.000000 

0.000000 

0.001415 

CO 

0.123704 

0.115658 

0.122957 

0.044893 

C02 

0.383614 

0.391660 

0.384361 

0.462446 

H 

0.000001 

0.000000 

0.000001 

0.000000 

H2 

0.054809 

0.062854 

0.055556 

0.128465 

H20 

0.437872 

0.429827 

0.437125 

0.362781 


MASS FRACTIONS 
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CH4 0.000000 0.000000 0.000000 0.000799 

CO 0.122237 0.114286 0.121499 0.044235 

C02 0.595582 0.608074 0.596742 0.715947 

H 0.000000 0.000000 0.000000 0.000000 

H2 0.003898 0.004470 0.003951 0.009110 

H20 0.278283 0.273170 0.277808 0.229909 


OADDITIONAL PRODUCTS 

ASSIGNED CONDITIONS 

WHICH WERE 

CONSIDERED BUT 

WHOSE MOLE 

FRACTIONS WERE LESS 

THAN .0000005 

FOR ALL 

C 

HCO 

CH 

H02 

C2 

C2H2 

C3 

C4 

JET-A(G) 

C12H26 

H202 

O 

OH 

02 

C (GR) 

H20(L) 

JP10(L) 

JET-A(L) 


C12H26(L) 


NOTE 

WEIGHT FRACTION OF FUEL IN TOTAL FUELS AND OF OXIDANT IN TOTAL OXIDANTS 
(SF) STANDS FOR (SHIFTING FROZEN) 

1 

0 FROZEN TRANSPORT PROPERTIES CALCULATED FROM EQUILIBRIUM CONCENTRATIONS 


CHAMBER 

THROAT 

EXIT 


MU 

(LBF-SEC/FT**2) 
1.09539496E-06 
1.02199533E-06 
6.93548600E-07 


K 

(LBF/SEC-DEG R) 
1.75734852E-02 
1.624 69596E-02 
1.09325694E-02 


VISCOSITY EXPONENT (OMEGA) FOR THE FORM MU=MUREF*(T/TREF) 1 


6.45086110E-01 
6.38496220E-01 
5.79323411E-01 
^OMEGA IS 0.76373 


MUREF FOR INPUT TO BLM= 

SPECIES CONSIDERED 
C 

C02 

H20 

02 


3.53317955E-05 LBM/(FT-SEC) 

IN TRANSPORT PROPERTIES CALCULATIONS 
CH CH4 CO 

C2 C2H2 H2 

H202 O OH 


ZONE = 1 

THEORETICAL ROCKET PERFORMANCE ASSUMING FROZEN COMPOSITION DURING EXPANSION 


0PC = 367.4 PSIA 


FROM AN ASSIGNED PRESSURE AND TEMPERATURE 

WT FRACTION ENTHALPY 


STATE TEMP DENSITY 

CHEMICAL FORMULA (SEE NOTE) CAL/MOL 

DEG K G/CC 

FUEL C 1.00000 H 1.94230 1.00000 -5430.000 

L 298.15 0.7730 

OXIDANT O 2.00000 1.00000 0.000 

G 298.15 0.0000 

0O/F=3.0000E+00 PERCENT FUEL=2.5000E+01 EQUIVALENCE RATIOl. 1343E+00 STOIC MIXTURE RATIO=3.4030E+00 
DENSITY=0.OOOOE+OO 


1 

CHAMBER 

THROAT 

EXIT 

EXIT 

PC/P 

1.0000 

1.7776 

1.0602 

43.153 

P, ATM 

25.00 

14.06 

23.58 

0.5793 

T, DEG K 

1500 

1359 

1485 

752 

H, CAL/G 

-1845.1 

-1902.7 

-1851.2 

-2133.5 

S, CAL/ (G) (K) 

2.0645 

2.0645 

2.0645 

2.0645 

G, CAL/GRAM 

-4941.9 

-4709.1 

-4917.5 

-3686.0 

U, CAL/GRAM 

-1950.2 

-1997.9 

-1955.3 

-2186.2 

DEN (G/LITER) 

5.76E+00 

3.57E+00 

5.48E+00 

2.66E-01 

M, MOL WT 

28.347 

28.347 

28.347 

28.347 

CP, CAL/(G)(K) 

0.4137 

0.4047 

0.4128 

0.3507 

GAMMA (S) 

1.2039 

1.2093 

1.2044 

1.2496 

SON VEL,M/SEC 

727.8 

694.4 

724.3 

525.0 

MACH NUMBER 

0.0000 

1.0000 

0.3124 

2.9602 

AE/AT 


1.0000 

2.0000 

6.0002 

CSTAR, FT/SEC 


3349 

3349 

3349 

CF VAC 


1.243 


1.661 

CF 


0.680 


1.523 

IVAC,LBF-S/LBM 


129.36 


172.94 

I, LBF-SEC/LBM 


70.81 


158.47 


MOLE FRACTIONS 


CO 0.123704 C02 0.383614 H 

0.054809 


0.000001 H2 


79 



H20 


0.437872 


MASS FRACTIONS 


CO 

0.278283 

0.122237 

C02 

0.595582 

H2 

0 . 

003898 

H20 


OADDITIONAL PRODUCTS WHICH 

ASSIGNED CONDITIONS 

WERE 

CONSIDERED BUT WHOSE 

MOLE FRACTIONS WERE 

LESS 

THAN .0000005 

FOR ALL 

C 

HCO 

CH 

H02 

C2 

C2H2 

C3 

C4 


JET-A(G) 

C12H26 

H202 

C12H26(L) 

0 

OH 

02 

C (GR) 

H20(L) 


JP10(L) 

JET-A(L) 


NOTE 

WEIGHT FRACTION OF FUEL IN TOTAL FUELS AND OF OXIDANT IN TOTAL OXIDANTS 

1 


****************************************************************************************************** 

****************** 

CALCULATE ODE AREA RATIO AND PRESSURE SCHEDULES FOR ZONE 1 

****************************************************************************************************** 

****************** 

*** EO f ENCOUNTERED IN READING ODE REACTANTS DATA *** 

*********CPU(SEC) = 0.0 ********** 
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